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FOREWORD 


oping and changing technologies and industrial products frequently 
d generation of materials that, if improperly dealt with, can threaten 
nvironment. The U. S. Environmental Protection Agency is charged 
| the nation's land, air, and water resources. Under a mandate of 

, the agency strives to formulate and implement actions leading to a 
| human activities aad the ability of natural systems to support and 
*ct the EPA to perform research to define our environmental prob- 
and search for solutions. 


ngineering Laboratory is responsible for planning, implementing, 
elopment, and demonstration programs to provide an authoritative, 
in support of the policies, programs, and regulations of the EPA 

er, wastewater, pesticides, toxic substances, solid and hazardous 

d activities. This publication is one of the products of that research 
ication link between the researcher and the user community. 


his document has application to the remediation of leaks from 
elines containing petroleum or other hazardous chemicals. 


E. Timothy Oppelt, Director 
Risk Reduction Engineering Laboratory 
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Experiments were conducted at the UST Test Apparatus Pipeline in which three acoustic 
sensors separated by a maximum distance of 38.1 m (125 ft) were used to monitor signals pro- 
duced by 11.4-, 5.7-, and 3.8-L/h (3.0-, 1.5-, and 1.0-gal/h) leaks in the wall of a 5-cm- (2-in.-) 
diameter pressurized petroleum pipeline. The line pressures and hole diameters used in the 
experiments ranged from 69 to 138 kPa (10 to 20 psi) and 0.4 to 0.7 mm (0.01 to .03 in.), respec- 
tively. Application of a leak location algorithm based upon the technique of coherence function 
analysis resulted in mean differences between predicted and actual leak locations of 
approximately 10 cm. The standard deviations of the location estimates were approximately 

30 cm. This is a significant improvement (i.e., smaller leaks over longer distances) over the 
cross-correlation-based techniques which are currently being used. 


Spectra computed from leak-on and leak-off time series indicate that the majority of acous- 
tic energy received in the far field of the leak is concentrated in a frequency band from 1 to 
4 kHz. The strength of the signal within this band was found to be proportional to the leak flow 
rate and line pressure. Energy propagation from leak to sensor was observed via three types of 
wave motion: longitudinal waves in the product, and longitudinal and transverse waves in the 
steel. The similarity between the measured wave speed and the nominal speed of sound in gaso- 
line suggests that longitudinal waves in the product dominate the spectrum of received acoustic 
energy. The effects of multiple-mode wave propagation and the reflection of acoustic signals 
within the pipeline were observed as non-random fluctuations in the measured phase difference 
between sensor pairs. 


Additional experiments with smaller holes and higher pressures (138 to 345 kPa [20 to 
50 psi]) are required to determine the smallest leaks that can be located over distances of several 
hundred feet. The current experiments indicate that improved phase-unwrapping algorithms 
and/or lower noise instrumentation are required to optimize system performance. 


This report was submitted in fulfillment of Contract No. 68-03-3409 by Vista Research, 
Inc., under the sponsorship of the U.S. Environmental Protection Agency. This report covers a 
period from 23 January 1991 to 31 October 1991, and work was completed as of 
30 September 1991. 
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Figure 1.1. Example of a passive-acoustic leak location system. ...................-.- ee 
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Figure 5.1. Diagram of the pressurized petroleum pipeline at the UST Test Apparatus. 
Pipe material is 5-cm- (2-in.-) diameter) steel; product is gasoline. Pressurized CO, 15 
used to generate 0- to 172-kPa (0- to 25-psi) static line pressure. Valves in connecting 


branches were closed during all experiments. 
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Figure 5.2. Diagram of the coupling between the acoustic transducer (CTI-30) an 
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Fisure6.1. Time series of acoustic leak signals generated by a 11.4-L (3-gal/h) | 
into a sand backfill. Line pressure is 103 kPa (15 psi) and hole diameter is 0.7 mm 
(0.03 in.). Sample rate is 10 kHz. A no-leak time series recorded by sensor B is sh 
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Figure 6.2. Time series of acoustic leak signals generated by a 1.9-L/h (0.5-gal/h) 
into a sand backfill. Line pressure is 345 kPa (5 psi) and hole diameter 15 0.3 mm ( 
in.). Sample rate is 10 kHz. A no-leak time series recorded by sensor B is shown f 
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1024-point data segments. Dashed lines indicate 95% and 99% levels of statistical 
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Figure 6.7. Coherence amplitude and coherence phase for acoustic leak signals 
bracketing a 1.9-L/h (0.5-gal/h) leak. Line pressure is 34 kPa (5 psi); sensor separat 
is 30 m (100 ft). Dashed lines indicate 95% and 99% levels of statistical significanc 
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included. The 
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Figure 7.1. Unwrapped coherence phase between 2.0 and 2.5 kHz for sensor pairs 
A-B, A-C, and B-C of Figure 5.1. Least-squares regression lines through actual data 
points are included. The flow rate is 11.4 L/h (3 gal/h). res. 
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delay (t) between time series recorded by sensors B and C. The time series were 
bandpass-filtered between 1.0 and 4.0 kHz prior to computing p,. The flow rate is 
11.4 L/h (3 gal/h). t. represents the predicted B-C time delay at V = 1000 m/s. .......... 
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between time series recorded by sensors B and C. The time series were bandpass- 


filtered between 2.0 and 2.5 kHz prior to computing p,. The flow rate is 11.4 L/h 


(3 gal/h). tc and ፒር; represent predicted time delays for primary and reflected 
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Figure 8.1. Time series of impulsive calibration signals recorded by sensors B and C 
of Figure 5.1. The estimated propagation speed (6250 m/s) is consistent with the 
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Figure 8.2. Unwrapped coherence phase between 2.7 and 3.0 kHz for sensor pair B-C 
of Figure 5.1 in which CO, is used as the product. The line pressure is 103 kPa 

(15 psi); the hole diameter is 0.7 mm. The estimated propagation speed (2400 m/s) is 
consistent with the nominal speed of transverse waves in steel. ሩሩ ssoseseseee. 


Figure 8.3. Unwrapped coherence phase between 2.1 and 2.4 kHz for sensor pair B-C 
in which the linear trend has been removed. The flow rate is 11.4 L/h ON. mee 


Figure 9.1. Unwrapped coherence phase between 1.5 and 4.5 kHz for sensor pairs 
A-B, A-C, and B-C. Solid lines indicate predicted coherence phase for linearly 
propagating plane waves based upon known leak location and propagation speed. 
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Figure 9.2. Unwrapped coherence phase between 1.5 and 4.5 kHz for sensor pairs 
A-B, A-C, and B-C. Solid lines indicate predicted coherence phase for linearly 
propagating plane waves based upon known leak location and propagation speed. 
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SECTION 1 
INTRODUCTION 


Underground pressurized pipelines are frequently used to transfer liquid products for many 
industrial applications. Some of these pipelines are associated with the underground storage 
tanks typically found at retail stations, and others with tanks at industrial storage facilities; they 
can contain petroleum products or a variety of other chemicals. There are many systems that can 
be used to detect leaks in underground pressurized pipelines. These leak detection systems are 
designed for use on pipelines that are typically 5 cm (2 in.) in diameter and generally 15.2 to 
61.0 m (50 to 200 ft) in length. The EPA regulations [1] require that the leak detection equip- 
ment used to test a pipeline monthly be capable of detecting leaks at least as small as 0.76 ፲ /ከ 
(0.2 gal/h) with a probability of detection (Pp) of 95% and a probability of false alarm (P,,) of 
3%. ፲፻ the equipment is used to test the line annually, 11 must be able to detect leaks as small 8. 
0.38 L/h (0.1 gal/h); in the regulations, this type of test is designated as a line tightness test. 


If a leak is found, remediation must follow. The first step in any remediation is to find the 
location of the leak. At the present time, there are two methods being used, but neither method is 
totally acceptable. The first method is to systematically uncover the line and perform a visual 
inspection for leaks. While this method works, it is time consuming, disruptive to operations, 
and costly. In addition, the line is subject to damage during the excavation process. The second 
method is to use a helium- or halogen-tracer technique, but both tracer techniques have opera- 
tional and accuracy problems. There is a need for a nondestructive method of leak location that 


15 accurate, relatively simple to use, and applicable to a wide variety of pipelines and pipeline 
products. 


One method of expediting the remediation process is the application of remote sensing 
techniques to the pipeline in order to accurately locate the leak. Passive-acoustic measurements, 
combined with advanced signal-processing methods, may provide a means by which to locate 
small leaks in limited-access pipeline delivery systems. The concept of using passive acoustics 
to determine the spatial location of leaks has been around for some time, but this approach has 
not been applied to underground petroleum pipelines. While it is known that a pressurized 
underground pipeline that is leaking emits an acoustic signal, the strength and characteristics of 
the signal associated with the leak are not well known. 








Acoustic systems have been successfully used to detect and locate leaks in nuclear reactors 
for many years [2]. By means of a cross-correlation analysis, 100- to 400-kHz acoustic sensors 
spaced at 5- to 10-m intervals can be used to detect leaks of about 3.8 L/min (1 gal/min) with an 
accuracy that is within 0.5 m. A similar approach has been tested for locating water leaks in 10- 
to 25-cm- (4- to 10-in.-) diameter underground district heating and cooling pipes [3]. Theoretical 
predictions based on Kupperman et al. [3] suggest that leaks of 7.6 L/min (2 gal/min) could be 
pinpointed to within several meters with sensors spaced at several hundred meters. Using moni- 
toring frequencies less than 25 kHz makes this wider spacing possible; frequencies between 1 


“and 5 kHz appear to give the best results. In the program of experiments conducted in 


Kupperman et al. [3], the sensors were externally mounted to the steel pipe. Interestingly, leaks 
that occurred in a steel pipe covered with insulation material (urethane and a rubber jacket) 
showed a higher level of signal intensity than leaks that occurred in an uncovered pipe. This 
suggests that if a pipeline is located underground rather than above ground, the surrounding 
backfill and soil may enhance the acoustic leak signal. The predictions were based upon models 
that were validated using a 25.9-m- (85-ft-) long, 7.6-cm- (3-in.-) diameter aboveground pipeline 
that contained water and in which leaks of 5.7 L/min (1.5 gal/min) were induced at a pressure of 
827 kPa (120 psi). Field tests on an operational line were also done. 


Figure 1.1 shows a simple representation of a passive-acoustic leak location system in 
which a pair of transducers bracketing a leak simultaneously sample the acoustic signal. The 
leak emits an acoustic wave, which propagates down a pressurized pipeline, and is received at 
the two transducers. ۸ microcomputer system acquires and analyzes the data needed to find the 
location of the leak. The data are converted from an analog signal to a digital signal with an A/D 
converter. These time series, recorded by spatially separated sensors, then serve as input to a 
leak location algorithm. 


Cross-correlation analysis works well provided that the signal is very strong or that the 
background noise is not excessive. When the acoustic signal is weak in relation to the level of 
background noise or has a finite frequency bandwidth, more sophisticated signal processing tech- 
niques are available. Advanced signal processing is required if any of the following objectives 
are to be achieved: (1) detection of leaks smaller than several gallons per hour, (2) a reduction in 
the number of false alarms and missed detections due to operational or ambient noise, and (3) an 
increase in ¿ne distance between sensors bracketing the leak. One such technique is coherence 
function analysis. 
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Figure 1.1. Example of a passive-acoustic leak location system. 


The application of coherence function analysis to signals measured by two or more trans- 
ducers is the means by which the source of ihe signal is best located. Coherence function analy- 
sis, which estimates the correspondence between two measurements as a function of frequency, 
is analogous to the squared correlation coefficient, but is a far more powerful tool in signal 
estimation and location. The coherence magnitude measures the strength of the correspondence, 
and the coherence phase measures the relative time delay. In contrast, the correlation coefficient 
15 a measure of correspondence that is the result of an integration over all frequencies. If the 
correspondence is frequency-dependent, or if the phase dependence of the correspondence is a 
nonlinear function of frequency, the correlation is degraded. By contrast, coherence is a direct 
measure of the complex frequency correspondence between two measurements, and, therefore, 
preserves the actual correspondence between the two measurements of the signal. 


In the last five to ten years there have been significant advances in commercially available 
acoustic sensors, in powerful computers that are both small and inexpensive, and in digital signal 
processing. This means that acoustic leak location can be made available in a portable package, 
a possibility that makes it an attractive and viable option. Acoustic systems are attractive from 
an operational standpoint because the test is short (a few minutes) and the sensors can be 
mounted directly on the outside of the pipeline. Acoustic systems have direct application to the 
15.2- to 61.0-m (50- to 200-ft) pipelines found at retail service stations because the sensors can 
be placed at each end of the line. Acoustic systems are also applicable to longer pipelines, but 
access at several points along the line must be provided. 


The objective of this work was to make an estimate of the accuracy of locating a leak in a 
pressurized petroleum pipeline, by means of passive acoustic sensors mounted on the outside 
wall of the pipeline, as a function of leak rate and distance between acoustic sensors. While 
there are regulatory standards for detection of leaks in underground pressurized pipelines, there 
are no standards for leak location. For rapid remediation, it would be highly desirable if the leak 
could be located within 1046 of the length of the pipeline in the case of a line longer than 30.5 m 
(100 ft), or within 3.0 m (10 ft) in the case of a line shorter than 30.5 m (100 ft). This limits the 
excavation to only a small fraction of the line. As will be shown below, theoretical estimates 
suggest that the accuracy of the proposed technique should be better than 25 cm. These theoreti- 
cal estimates, however, assume that the leak signal is large compared to the noise. 


It would be desirable if the acoustic system could locate leaks as small as the detection 
standards in EPA”s regulations (i.e., either 0.38 or 0.76 L/h [0.1 or 0.2 gal/h]). However, the 
system can be useful nonetheless, even if it does not meet these standards, since many of the 
leaxs actually detected are much larger. One to two drops of liquid per second produces leaks as 
large as 0.38 L/h (0.1 gal/h); such a leak might occur if, for example, the threaded connection 
between two pipes or between a pipe and another appurtenance is not tight. Any type of small 
hole in the pipeline system will produce a large leak. A 0.08-cm- (1/32-in.-) diameter hole in a 
pipeline pressurized to 138 kPa (20 psi) would, for example, result in a flow rate of approxi- 
mately 19 L/h (5.0 gal/h). Leaks of this magnitude should be large enough to locate with an 
acoustic system providing that the length of the line is not too long. 


The results of this work are described below and are summarized in a peer-reviewed jour- 
nal art.cle being published by the American Society for Testing and Materials [4]. 


SECTION 2 
CONCLUSIONS 


Passive acoustic measurements, combined with advanced signal processing techniques 
based on coherence analysis, offer a promising method for the location of small leaks in pressur- 
ized petroleum pipelines found at retail service stations and industrial storage facilities. While 
the results presented in this work represent a significant improvement over previous pipeline leak 
lecation efforts, additional research and development are required in order to optimize system 


performance. Location of leaks of several tenths of a gallon per hour over distances of several 
hundred feet should ultimately be possible. 


Experiments were conducted on a pipeline at the UST Test Apparatus in which three 
acoustic sensors separated by a maximum distance of 38 m (125 ft) were used to monitor signals 
produced by 3.8-, 5.7-, and 11.4-L/h (1.0-, 1.5-, and 3.0-gal/h) gasoline leaks. These flow rates 
were generated through drilled holes 0.4 to 0.7 mm (0.01 to 0.03 in.) in diameter. The three- 
transducer system enabled the propagation speed of acoustic waves to be measured for particular 
combinations of product, pipeline geometry, and analysis frequency band. Data recorded at the 
higher flow rates (5.7 and 11.4 L/h [1.5 and 3.0 gal/h]) correspond to full line pressure (103 to 
138 kPa [15 to 20 psi]), while data recorded at the lower flow rate (3.8 L/h [1.0 gal/h]) were 
obtained under partial line pressure (69 kPa [10 psi]) due to the limitation imposed by the mini- 
mum available hole diameter (0.4 mm [0.01 in.]). Application of a leak location algorithm based 
upon the technique of coherence function analysis resulted in mean differences between 
predicted and actual leak locations of 8.7 cm (11.4 L/h [3.0 gal/h]), 3.6 cm (5.7 L/h [1.5 gal/h]), 
and -11.6 cm (3.8 L/h [1.0 gal/h]). Standard deviations of the location estimates were 26.1 cm 
(11.4 L/h [3.0 gal/h]), 26.3 cm (5.7 L/h [1.5 gal/h]), and 39.1 cm (3.8 L/h [1.0 gal/h]). The mean 
propagation speed was 915 m/s with a standard deviation of 146 m/s. 


Data recorded in the presence cf a 1.9-L/h (0.5-gal/h) leak were obtained as part of an 
investigation of signal strength as a function of line pressure for a fixed-diameter hole (0.4 mm 
[0.01 in.]). The 1.9-L/h (0.5-gal/h) leak produced a detectable signal; however, because of the 
reduced line pressure, the algorithm, as applied, yielded no location estimates. 


Spectra computed from leak-on and leak-off time series indicate that the majority of acous- 
tic energy received in the far field of the leak is concentrated in a frequency band from 1 to 
4 kHz. The strength of the acoustic signal within this band was observed to be proportional to 
the leak flow rate and line pressure, as expected. Energy propagation from leak to sensor was 
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separations. Asa conseguence, the measured phase cannot be accurately unwr: 
within frequency bands where y?(f) is high; thus, a differential form of Eq. (4.3 
to relate sensor separation, propagation speed, and coherence phase: 


leak rates and at one or more separation distances of the sensors bracketing the leak. The algo- 
rithms proposed for leak location and the estimates of the uncertainty in applying these algo- 
rithms are described below. 


4.1 Location Algorithm 





f wave motion: longitudinal waves in the product, and both transverse 
the steel. Isolation of each of these propagation modes was achieved 
and CO, as the product fluids, and through the generation of impul- 
hough each of these propagation modes is believed to contribute to the 
¡gitudinal wave motion in the product was clearly the dominant 
d-filled pipelines. The effects of multiple-mode wave propagation 

DC signals within the pipeline were observed as non-random fluctua- 

> difference between sensor pairs. 














SECTION 4 
LOCATION OF A CONTINUOUS LEAK SIGNAL 


SECTION 3 
RECOMMENDATIONS 















































































The primary function of the location algorithm is to estimate the time delay between 
acoustic leak signals received by a pair of sensors. The time delay, which is estimated from 
phase measurements made in the frequency domain for continuous leak signals, can be used to 
estimate the source location (for sensors bracketing the leak) or the propagation speed of the 


The full capability of the location algorithm was not evaluated in these tests. The smallest 
hole used to generate a leak in the experiments was 0.4 mm (0.01 in.). Ata line pressure of 
138 kPa (20 psi) this resulted in a leak rate of 3.8 L/h (1.0 gal/h). It is recommended that addi- 
tional experiments be performed with smaller holes at higher line pressures (138 to 345 kPa [20 
to 50 psi]) to determine the minimum leak rate that can be reliably located. The current work 
indicates that further improvement can be realized through the application of better phase- 
unwrapping algorithms and better instrumentation. A better understanding of the underlying 
physics of pipeline acoustics, including the propagation modes and source mechanisms of the 
acoustic leak signal, will help optimize the algorithms and the hardware. It is recommended that 


Consider two measurements of the signal, m,(t) and m,(t), where each represents the 
sum of a desired signal, s(t), and a contaminating noise component, n(t). The signal, s(t), 
would be the acoustic signal emanating from the leaking pipeline, and the contaminating 











in which, for the sake of simplicity, it is assumed that the medium is nondispers 







The three-sensor approach illustrated in Figure 4.1 was used to locate leal 
underground pipeline. This configuration is similar to the one shown in Figure 
that sensor pair B-C is used to measure the in situ wave speed, while sensor paii 
are used to estimate the leak iocation. Knowledge of the wave speed improves ! 
of the leak location estimate. However, the wave speed associated with a partic 
and pipeline geometry is usually unknown. Therefore, it is useful to make an e 






noise component, n(t), could be ambient noise in the measurement environment that is uncor- 
related at the separated acoustic sensors. The coherence function, Y (f), is the normalized 
cross spectrum of the two measurements, 


acoustic waves (for non-bracketing sensor pairs). Two criteria must be satisfied in order that 





n requires the identification of frequency bands within which a high 
itained between acoustic signals propagated along different paths from 
function analysis provides the best means of gauging this similarity, 




























accurate location estimates result from the application of the location algorithm: (1) the received 
signals must originate primarily at a single, localized source and propagate as plane waves along 

































(or within) the pipeline, and (2) the received signals must maintain a reasonable degree of simi- A= M,(f)M,(f) (41) 
VIM AP VIAL 


information concerning the leak location from ambient or system larity over the maximum sensor separation. If criterion (1) is satisfied, the difference in phase 





noise ratio (SNR) was observed to be generally high within the entire the following work be performed to extend the technology: between received waves of a given frequency is simply related to the time delay between signals estimate of the wave speed. 
|, continuous regions of high coherence appropriate for source loca- 07 : : i that arrive at the different sensor locations. The accuracy with which the time delays can be where upper-case letters denote the Fourier transform of the respective quantities and the A LEAK B e 
00 Hz in width. Several data sets recorded in the presence of the . et: > መመ z b em سم‎ ` capable of processing the coherence phase over measured is related to criterion (2). The similarity between signals emitted from a localized overbar denotes the ensemble average. As noted above, the coherence is complex and the 4 - | 


phase, $(f), measures the relative time delay between the two signals at a given frequency f. 
The coherence function ranges in magnitude from 0 (signals completely uncorrelated) to 1 
(signals completely correlated). Values of Y(f) exceeding 95% of the noise fluctuations are 
usually taken as indicating a reliable phase measurement. 


ibited high coherence over a 2-kHz bandwidth. Location estimates 
‘oss-correlation technique showed that without the detailed knowl- 
vided by the coherence function, cross-correlation analysis cannot 


source and received at separate locations is determined by the signal strength relative to ambient 
noise (i.e., the signal-to-noise ratio) and the difference in propagation path between the source 
and each sensor. Due to the complex manner in which the acoustic leak signal is produced (tur- 


* characterization of the wave propagation modes excited by the acoustic leak signal 
and the degree to which each mode enhances or degrades the leak location estimate 


* reduction of system noise through transducers specifically designed for low- 
frequency, high-sensitivity applications, and through the use of low-noise, audio- 
range preamplifiers 
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ptable accuracy. The observed correspondence between measured 


vithin the 1- to 4-kHz frequency band demonstrates the need to 
d location algorithm such that a greater fraction of the info:mation 
ignals may be processed. 


bulent flow and cavitation) and the many variations in the propagation medium (valves, 















branches, reflective ends), the degree of signal similarity is not uniform over a broad range of If the acoustic leak signal is approximated as a collection of propagating acoustic plane 


* automation of the data acquisition system and signal processing algorithm, and evalu- idit diy dii ከለ ےش‎ de i lati 
waves that obey the simple linear dispersion relation 


Figure 4.1. Three-sensor approach to acoustic location of leaks. 
ation of system performance on a variety of actual pipelines 


Application of Eq. (4.4) to sensors A and B, which bracket the leak, yields 
relationship between measured phase, wave speed, and leak location: 


frequencies. Though the signal-to-noise ratio provides a reasonable estimate of the frequency 
band for which accurate leak locations may be obtained, a more sensitive measure of Signal simi- 
larity is required for the location of small leaks (e.g., 11.4 L/h [3 gal/h] or less). 
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ide a generally quiet ambient environment in which to perform 




















۱٥١ the SNR for a given leak largely determines the ability of a pas- The proposed detection algorithm is designed to locate leaks with a high signal-to-noise where K is the wavenumber and V is the propagation speed, the differential separation = Kas لچ‎ Ada 
te the leak, the system noise level should be determined by ambient ratio'. As the leak becomes smaller or the line becomes longer, the SNR will eventually become between two sensors, Ar, and the frequency-dependent phase, $(f), are simply related by — di 

lectronic noise. The combination of sensors (CTI-305) and small enough that routine location estimates cannot be made. The accuracy of the location esti- 00) = 2x Ar ) (4.3) Xo V ፀቆህ 
5660-Cs) used in this work was incapable of resolving the low levels mate depends on the estimate of the sound speed, the resolution of the data acquisition system, V ۵ 47 ۴ 


sociated with the pipeline at the UST Test Apparatus. Improved sys- 
ained through the use of transducers with greater sensitivity in the 
) kHz) and low-noise preamplifiers. 


the strength of the signal compared to the noise, and the signal processing algorithm. The 
experimental uncertainty in locating a leak can be estimated from the standard deviation of the 
location estimates made from many realizations of the measurement of the leak for one op more 


Through the use of coherence function analysis, it is possible to isolate portions of the acous- 
tic spectrum within which the linear dispersion relation is obeyed. The measured phase shift, 
$(f), within these frequency bands can then be used to estimate either the propagation speed 
of acoustic waves or the differential sensor separation. Because the coherence phase is con- 
fined to the range -180° < ቀ < 180°, the measured phase generally differs from the actual 
phase by an unknown factor of 360°, except at very low frequencies and/or very small sensor 


where L denotes the location of the leak. The wave speed 15 estimated from the 1 
phase between sensors B and C: 
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1 The signal-to-noise ratio is estimated from the ratio of the received power at the output of the measurement system 
over a specified frequency band while a signal is present, divided by the output power measured at the receiver over 
the same frequency band while a signal is not present. 
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»pped except 4.2 Location Errors o{X,,} = <V ا‎ dé | በ +ጽጓ” | (412) Figure Pan shows the uncertainty in phase determined from Eq. (4.9) as 8 function of 1 
) must be used — 47 df the number of independent data segments, n, used to compute the coherence and for ¥;,,(f) = 
The coherence phase required for the estimation of the leak location and wave speed is 0.2, 0.4, 0.6, and 0.8. When n = 15 (which is the number of independent data segments used o 1 
most accurately measured (i.e., has a low probability of being a noise fluctuation) when the in the leak location estimates in this paper), and Y ¿(f) = 0.5 (which is typical of the values of m 28 - 
(4.4) coherence amplitude, (f), is significantly above the noise. The value of the coherence ` ESOS ار‎ coherence for many of the leak location estimates), 016, 5(/)) = 10.5” 9) 
amplitude associated with a level of statistical significance, P, is given by E? X2. 17177 0 i RES i : ۱ c 
Figures 4.3 through 4.5 show the uncertainty in the location estimate, G{X,,}, when the G 
ሦ=1--(1--ዖ)', B= 05m - D-1 (4.8) The uncertainty, o{X,, }, in the estimate of the location of the leak, يی‎ is affected by distances between sensors A and B and sensors B and C are nominally 38.1 m (125 ft) and 0 
— the three-sensor measurement geometry and the location of the leak. The term K accounts 7.6 m (25 ft), respectively, and the distance between the leak and sensor A is 3.8 m (12.5 ft), < 
1.1, except where n is the i;umber of independent data segments used to compute Y (f). If P is 0.95 and n for this effect; clearly, the largest uncertainty occurs when K is largest. Eqs. (4.12) and 9.7 m (31.75 ft), and 19.0 m (62.5 ft). The curves are derived in terms of the number, n, of a 
ዓል መሎ is 15, then for any value ' (f) greater than 0.35, there is less than a 0.05 (5%) chance that (f) (4.13) suggest that the minimum uncertainty in XAr occurs when the leak is midway between independent data segments used to estimate the coherence and ,(f). The estimates assume 
the accuracy was produced by a noise fluctuation. A and B; in fact, when X,, = 0.5 مبلا‎ the uncertainty in X,, is independent of the uncer- that d@/df is computed from N = 26 points over the frequency band 2119 to 2373 Hz; the 
— An estimate of the one-standard-deviation uncertainty in the phase measurements, tainty, o1d6gc/df), in the measurement of the propagation velocity. The maximum uncer- 254-Hz band represents the larger bands for which the location estimates were made. In 
(perimental 6{%,»(/)}, derived from coherence measurements between any two sensors (e.g., A and B) is tainty in X, occurs when the leak is very near sensor A (or B), or when the separation many instances a frequency band closer to 100 Hz was used. The conditions used to generate NUMBER OF INDEPENDENT DATA SEGMENTS 
given by distance, Kpc, used to estimate the propagation velocity is very small. Figure 4.6 are identical to those used to generate Figure 4.4, except that N = 5 points. Figure 4.2. 0:19,,(/)] for sensor separation distance asa function of n between 5 and 50 and values of وو‎ 
]1- 7 The estimates of uncertainty, G{X,, ን, can be written in terms of the standard deviation | ge — of the leak tê and sensor separation is best made — * value of equal to 0.8, 0.6, 0.4, and 0.2 estimated for any frequency with Eq. (4.9). 
— Obs} = — — (4.9) Y (f) is significantly above the noise. Table 4.1 shows the one standard deviation in phase, 
Vs Van of the phase, 0{$}, using the equation sensor separation, and leak location as a function of n for the three geometries described 
by Bendat and Piersol [5]. Eq. (4.9) indicates that the random error associated with phase s zi .. _ 7 (414) above for values of Get) at the 95% and 99% levels of statistical significance. As n 
measurement is controlled by n and by the magnitude of ጄኔ). dí) NEG increases, the value of 45 (f) required to maintain the same level of statistical significance E 
The 06-68118) ተከ direi Gib: دا‎ eibi di مع‎ af a هده‎ | | i e | decreases. Obviously, more accurate estimates of the location of the leak can be made if the T 
y, ALS» E e where N is equal to the number of independent data points in the frequency band, Af, used to values of PaM or n are greater than the values in Table 4.1. The results in Figures 4.3 5 
co with get (0 መዋ ላ, when sensors A and 8 A - leak and the mean propagation — the rate of change of the phase. In Am" there is some correlation between data through 4.5, as well as in Table 4.1, suggest that the largest uncertainty in leak location is less 7 
velocity, <V>, is measured between sensors B and C, is given by points, and so the number of degrees of freedom is usually less than the number of data than 20 cm (7.88 in.). Z 
j ۱ lé dt, ý ۱ 22 points used to derive 10/1۸ Estimates of 0{X,,} are made in the next section for different 7 
(4.5) o{X,,} = um s u) + A s «ቁረ | | (4.10) sensor geometries, levels of statistical significance, and frequency bands. 
| M \ " ri > | 4 j The most accurate estimates of leak location are made when (1) the number of indepen- 
(4.6) Because the average rate of change of the phase with respect to frequency is directly related dent data sets, n, used to IER wa ^ nce and phase, the number of degree +e freedom, N, 
to physical distances, it can be shown that and the frequency bandwidth used to estimate do/df are large, and (2) the leak is 
approximately at the mid-point between the two bracketing sensors. 
neasured NUMBER OF INDEPENDENT DATA SEGMENTS 
4.11) 90 0 یر سس ھی و‎ gagap eme رت سے و عو میڈ‎ 
(47) A theoretical estimate of the standard deviation of the location estimate, O(X,, ), can be ARUM candhana: Ñ 5 سو‎ Ett an^ m (25 ft), ጁኒ = 3.8 m (12.5 ft) (i.e., 10% of Nal, Af = 


Assuming that o(d(/df) = o(d$,,/df) = G{dQgc/df}, and using Eq. (4.11), Eq. (4.10) 
reduces to 


made with Eqs. (4.9), (4.12) and (4.14). A field estimate of O(X,,) can be made, from repet- 
itive measurements, with Eq. (4.10). 
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NUMBER OF INDEPENDENT DATA SEGMENTS 
Figure 4.4. 0{X,,} for sensor separation distance between two sensors, A and B, as a function of n between 5 


and 50 and values of 2, (f) equal to 0.8, 0.6, 0.4, and 0.2 estimated with Eas. (4.9), (4.12), and (4.14) for the 


following conditions: X,» = 38.1 m (125 ft), Ka = 7.6 m (25 ft), X,, = 9.7 m (31.75 ft) 72 e 
2373 - 2119 = 254 Hz, N = 26, and «V» = 1000 m/s O D» Xu, = 9.7 m (31.75 À) (i.e., 25% of Kan), Af 
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NUMBER OF INDEPENDENT DATA SEGMENTS 
Figure 4.5. 5{X,,} foi sensor separation distance between two sensors, A and 8, as a function of በ between 5 


and 50 and values of R.M equal to 0.8, 0.6, 0.4, and 0.2 estimated with Eqs. (4.9), (4.12), and (4.14) for the 


following conditions: X,, = 38.1 m (125 ft), Xy. = 7.6 m (25 ft), X,, = 190 m (62.5 ft) በ = 
2373 - 2119 = 254 Hz, N = 26, and cV» = 1009 m/s dr Ku m (62.5 ft) (i.e., 50% of አፌ), Af 








DISTANCE ERROR - cm 
5 58 6ይ 














Li 
5 15 25 35 45 


NUMBER OF INDEPENDENT DATA SEGMENTS 
Figure 4.6. 0(X,,) for sensor separation distance between two sensors, A and B, as a function of n between 5 


and 50 and values of 4, (/) equal to 0.8, 0.6, 0.4, and 0.2 estimated with Eqs. (4.9), (4.12), and (4.14) for the 


following conditions: X,» = 38.1 m (125 ft), Xgc = 7.6 m (25 ft), ጂነ 29.7 m (31.75 ft) (i.e., 25% of Kap), Af = 
2373 - 2119 = 254 Hz, N = 5, and «V» = 1,000 m/s. 


4.4 Accuracy 


The accuracy of the acoustic measurement system for locating a leak along the pipeline 
is defined as the mean difference between the actual and estimated location of the leak. The 
accuracy is affected by any biases that may result from the phase measurement estimates, by 
the uncertainty in the estimate of the average propagation speed (as affected by an error in 
the location of the sensors with respect to each other), and by differences in the wave speed 
along different propagation paths from the leak to the sensors. 


A bias in the phase measurements can occur for a variety of reasons: biases in the 
power spectral and cross spectral estimates used to compute phase, ambient noise that is 
highly directional, and other signals that are correlated with the input signal. These may 
include errors due to multiple modes of propagation that occur in the liquid product and the 
pipe wall and errors due to reflections of the signal within the pipe due to boundaries, joints, 
couplings, elbows, etc. Accurate theoretical estimates of bias errors are difficult to make 
and require information about the signal and noise field; as a consequence, none will be 
attempted here. Itis aniicipated, however, that the bias errors in the phase measurements will 
be small compared to the other two sources of error mentioned above. 


























Table 4.1. Estimates of the Total Uncertainty in Measuring the Separativa between Sensors and the Location of 
a Leak as a Function of the Number of Incoherent Averages and Y. ¿(f) Computed for L = 95 and 99% 
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"Sensor separation error estimated between 2119 and 2373 Hz for V = 1,000 m/s and N = 26 assuming that the 
mp | is measured and the nominal sensor separations are X,» = 38.1 m (125 ft) and Kass 7.6m (25 ft), and 
Kan = 9.7 m (31.75 ft). 

- Leak location error estimated between 2119 and 2373 Hz for V = 1,000 m/s and N = 26 assuming that the 
velocity is measured and the nominal sensor separations are X,» = 38.1 m (125 ft) and X,- = 7.6 m (25 ft), and 
X, = 3.8 m (12.5 fi). 

"Sensor separation error estimated between 2119 and 2373 Hz for V = 1,000 m/s and N = 26 assuming that the 
velocity is measured and the nominal sensor separations are X,» = 38.1 m (125 ft) and Xpc = 7.6 m (25 ft), and 
Kar = 9.7 m (31.75 ft). 

“Leak location error estimated between 2119 and 2373 Hz for Y = 1,000 m/s and N = 26 assuming that the 


velocity is measured and the nominal sensor separations are X,, = 38.1 m (125 ft) and Xy. = 7.6 m (25 ft), and 
X aL = 19.0 m (62.5 ft). 
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Ultimately, the accuracy of a leak location estimate depends on the speed of sound 
used to covert the time-of-arrival or phase measurements to distance. Any errors in the esti- 
mate of the wave propagation speed will translate directly into distance errors. Sound speed 
estimates should be measured experimentally to minimize these errors. 


4.5 Performance of the Leak Location System 


A leak detection system requires a high probability of detection (P,) to be effective. 
The EPA regulation requires that a detection system have a P, of at least 95%. On the other 
hand, a leak location system does not require a high Pp, but does require a low probability of 
false alarm (P,,), 5% or less, to be effective. There is no recommended Pr, or P, for a loca- 
tion system. However, a Pb of 50% or less should be sufficient because it is presumed that 
te leak location system is being used only after the existence of a leak has been confirmed. 


The PFA and Pp can be estimated directly from the SNR for normally distributed noise. These 
estimates can be made with the standard curves shown in Figure 4.7; these curves were 
reproduced directly from Urick [6] and can be found in many textbooks. 
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* 4.7. Estimates of P, and P,, as a function of SNR = d developed for normally distributed noise. (Urick 


Table 4.2 presents, as a function of SNR, estimates of the Pra for a Pp of 50% and Pp 
for a Pp, of 5%. For most detection applications, an SNR over 9 will more than suffice. The 
ሾለ and Pp can be further improved by averaging the results of multiple tests. 


Table 4.2. Estimates of P,, and a P, of 50% as a Function of SNR 





P,, for P, for 
SNR Po = 50% Pe, = 5% 


(%) 





50 
16 
ህ 
2.5 
0.15 
0.008 
< 0.0001 





The SNR can ኩይ computed from the coherence estimates made over the signal band 0፻ 


interest for any two sensors bracketing a leak in the pipeline [7]. The relationship between 
Y (f) and the SNR is given by 


SNR = -፻0) (4.15) 
1 -Y(f 

Eq. (4.8) can be used to estimate the probability of a false alarm as a function of SNR 
and the number of independent data segments used to compute the coherence. Thus, with 
Figure 4.7 (or Table 4.2) and Eqs. (4.8), and (4.15), (1) an estimate of the Pg, and Pp can be 
made for a leak location test based on the value of n used in the analysis and the value of (f) 
measured, or (2) the minimum specifications of n and Y (f) required to achieve a given Pe, 
and Pb can be estimated. The latter can be done by | 


e selecting the desired Pr, and P,, and computing SNR from Figure 4.7, 


* solving Eq. (4.15) for Y (f) and computing Y(f) for the SNR required to achieve the 
. desired Pr, and Pp, and 


computing n from Eq. (4.8) for Y (f) and P,,. 


The following example provides an estimate of the minimum number of independent data 
segments, n, and the minimum coherence, ÿ(f), required to achieve a given performance. If 
a Pra of 5% and a Py of 50% are desired, the SNR estimated from Figure 4.7 must be approx- 
imately 3 (6 dB). The value of (f) estimated from Eq. (4.15) is 0.75, and the value of n 


estimated from Eq. (4.8) is 5. If n is increased and the value of Y (f) remains at 0.75, the Pra 
will be less than 5% and the P, will greater than 50%. 


The Pp and رم‎ can be estimated for a given measurement by 


* computing the Pr, for a given n and the measured Y (f) by means of Eq. (4.8), 


* computing SNR for the measured Y(f) by means of Eq. (4.15), and 
* computing Pb for the given P- and SNR by means of Figure 4.7. 
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SECTION 5 
EXPERIMENT DESIGN 


lucted on the pressurized 5-cm- (2-in.-) diameter steel pipeline 
ere done in accordance with the quality assurance project plan 
the UST Test Apparatus is shown in Figure 5.1. Access ports 
1sducers to the pipeline were located at intervals of approxi- 
etween the pipeline and an acoustic transducer is shown in 

vn in Figure 5.1 were used during all experiments. The 

were CTI-30 resonant sensors. Though the CTI-30 is designed 
pplications, its sensitivity at low frequencies (1 to 5 kHz) is 
istic leak signals in pipelines. The acoustic signals were ampli- 
leans of battery-operated Panametrics 5660-C preamplifiers 
amplifying filters. A Western Graphtec TDA-3500 transient 
icoustic wave forms at a sampling rate of 10 kHz. Data were 
>AQ-386 portable computer. A diagram of the data acquisition 


of the apparatus used to generate a leak in the pipeline. The 

d by the static line pressure (0 to 172 kPa [0 to 25 psi]) and the 
which the product was allowed to leak. Leak apertures between 
in.) were introduced into the pipeline via carburetor jets in 
11በ፪ small-diameter holes through the steel wall of the pipe- 
rated during the experiments was between 1.9 and 18.9 L/h 
naterials used in the experiments were fine-grain sand and 


surements (calibration, leak-on, and leak-off) were made for 

, hole diameter, and backfill material. The calibration signal 

] lead on the pipe surface near the location of the simulated 
this impulsive signal at the three transducer locations were 
data acauisition system were operating properly. After the 

ly eight leak-on measurements, each 1.7 s in duration, were 
ight leak-on measurements were bracketed by a pair of record- 
tions. 
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Figure 5.4. Apparatus used to generate simulated pipeline leak. Backfill materials are fine-grain sand and coarse 


Figure 5.2. Diagram of the coupling between the acoustic transducer (CTI-30) and the steel pipeline. 
E" Leak apertures between 0.3 and 0.8 mm (0.01 to 0.03 in.) were introduced into the pipeline via carburetor 
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Figure 5.1. Diagram of the pressurized petroleum pipeline at the UST Test Apparatus. Pipe material is 5-cm- 
(2-in.-) diameter steel; product 15 gasoline. Pressurized CO, 15 used to generate 0- to 172-kPa (0- to 25-psi) static 
line pressure. Valves in connecting branches were closed during all experiments. 
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Figure 5.3. Diagram of the data acquisition system used in the experiments. 









SECTION 6 
DATA 


The raw data consist of time series of acoustic leak signals and ambient noise san 
simultaneously by three sensors. The first step toward deriving and applying a leak loc. 
algorithm to the raw data 15 to view the data in three forms: (1) time series, (2) power s 
density, and (3) complex coherence. Viewed in the time domain, the leak-on and leak-« 
(i.e., time series) provide convincing evidence that an acoustic leak signal exists and is : 
over the dimensions of the pipeline. In addition, the time series reveal the continuous c. 
of the leak signal, as compared to the impulsive signal generated by the pencil-lead cali 
signal. However, the time series alone offer no clues as to the location of the leak or the 
processing required to perform a source location estimate. The distribution, with respec 
quency, of acoustic energy emitted by the leak and the way in which this energy is prop 
from source to sensor is revealed by viewing the data in the frequency domain (i.e., pou 
spectra and complex coherence). 


Time series of acoustic leak signals generated by a 11.4-L/h (3-gal/h) gasoline lea 
sand backfili are shown in Figure 6.1; a time series recorded under no-leak conditions b 
the sensors is shown for reference. Aside from an anti-alias filter applied to the analog : 
prior to digitization, the time series shown here represent unfiltered data. The distance | 
the simulated leak and sensors A and B is approximately 15 m (50 ft); sensor C is locate 
approximately 23 m (75 ft) from the leak. The line pressure used in this experiment was 
(15 psi) and the hole diameter was 0.7 mm (0.03 in). Two important observations are m 
regarding the time series of Figure 6.1: (1) the leak is clearly detectable from the differe 
between the leak-on and lea”-off measurements, and (2) the relative arrival time of the 1 
nal at the different sensor locations cannot be obtained through inspection of the time se 
continuous nature of the acoustic leak signal requires that some type of signal processin; 
applied to the leak signal time series in order that the relative arrival times, and hence th 
tion of the leak, can be estimated. | 


Figure 6.2 shows time series of acoustic leak signals generated by a 1.9-L/h (0.5-g 
leak into a sand backfill along with a no-leak time series for reference. The line pressure 
this experiment was 34 kPa (5 psi) and the hole diameter was 0.3 mm (0.014 in.); the ser 
geometry is the same for all experiments. At the reduced flow rate of 1.9 L/h (0.5 gal/h) 
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acoustic leak signal cannot be clearly identified in the raw time series. In order to detect the : Ge., 1.9 L/h [0.5 gaV/h]). The SNR at each flow rate was obtained by dividing the power | پر ې‎ —QÜ ር d * TM "^D c 9 den e e = ووي‎ rA sat 
presence of small leaks, the data collected in the time domain must be transformed into the fre- diy spectral density computed with the leak present by a similar spectrum computed with no leak e = — d وپ = ` و‎ v) id a ۸ 
— dene. EE چی شف‎ present. The power spectra for each of the three individual sensors, computed using 31 over- | signal similarity is maintained is narrowed considerably. 
| lapping, 1024-point FFT segments, were averaged together prior to computing the SNR. The If the coherence amplitude is statistically significant for each Fourier component within 
ipled P LLLA MEL NN Md Alo Il A Al illl DURANT | time series used were 1.7 s in duration and were Sampled at a frequency of 10 kHz. The SNR a given frequency band, a simple phase-unwrapping procedure may be applied to the coher- 
0 i spectra show that the energy associated with the acoustic leak signal is not equally distributed | ence phase. Figure 6.8 shows the unwrapped phase shift as a function of frequency between 
"— ve LEAK OFF ae «٠ | | over the 1- to 5000-Hz sampling bandwidth, but is instead concentrated within a relatively 3.8 and 4.0 kHz for sensors 8 and C of Figure 5.1; the flow rate is 11.4 L/h (3 gal/h). 
ff data m D | T. | | | | narrow 1- to 4-kHz frequency band. The frequency domain representation of acoustic data 
— | D | offers a means by which the location algorithm can separate useful information concerning 
haracter y AAA LU LAS ABD LA auk da AE LI PT f AN LA Un ዜሰ፤ ሊዘ በቢ ITET, dau | the leak from unwanted noise. 
bration EE Figure 6.4 shows the SNR computed from time series recorded in the presence of an 
> types of 11.4-L/h (3-gal/h) leak with coarse gravel as the backfill material. A comparison of 
t to fre- Figures 6.4 and 6.3(a) indicates that the acoustic signal produced by the turbulent flow of 
agated product into a backfill of fine-grain sand represents a slight enhancement over that produced 
‚er when the backfill consists of gravel. A similar relationship between the backfill material and 
| the strength of the acoustic leak signal has been reported for studies of simulated leaks in the 
TW di | | | floor of aboveground storage tanks [9-11]. 
y one of An estimate of the strength of the acoustic signal as a function of line pressure for a 
signals fixed-diameter hole is shown in Figure 6.5. The line pressures used in this experiment were 
ግመ 172 kPa (25.0 psi), 128 kPa (18.5 psi), 83 kPa (12.0 psi), and 41 kPa (6.0 psi); the hole diam- 
d eter was 0.5 mm (0.020 in.). The average SNR computed over the 1- to 4-kHz frequency 
: 103 kPa h band 15 used as the measure of acoustic leak signal strength. As expected, the leak signal 
ade mu | | 1 |... ۸ : ۱ du. 100.0 150.0 increases in strength as the line pressure is increased. 
nce | — bii HL Lu c LO ab Lud il | | | TIME (mS) 
1 | ۱ ۱ i ۱ 6.2 Coherence and Phase Measurements 
eak sig- Figure 6.2. Time series of acoustic leak signals generated by a 1.9-L/h (0.5-gal/h) leak into a sand backfill. Line 
ries. The 500 100 0 50. 0 | | E Rte no 0 0.3 mm (0.01 in). Sample rate is 10 kHz. ላ no-leak time series EN Figure 6.6 shows the coherence amplitude and coherence phase as a function of fre- 
: be TIME (mS) d 6.1 Signal Strength | quency for — we signals received by sensors bracketing a 5.7-L/h (1.5-gal/h) leak. 
— Figure 6.1. Time series of acoustic leak signals generated by a 11.4-L/h (3-gal/h) leak into a sand backfill. Line | The strength of the acoustic signal led pipeline is primari ےہ دود د وو سیا — ې وور ور د‎ en pete ووو‎ 
| የመጽ is 10) RP (15 pel) cad ኩይ ak signala generated by à 11. San C MA) dk into a sand backfill. Lino ! m gnal produced by a leak in a — pipeline is —— of 15 overlapping, 1024-point segments, each individually detrended and weighted with a 
ሰ recorded by sensor 8 is shown for reference. | a 202 of — rate, but - — — by ms surrounding backfill material. cosine bell prior to Fourier transforming. Statistically significant coherence (as indicated by 
۱ Estimates of the signal-to-noise ratio for pipeline leaks into a sand backfill at flow rates of the 95% confidence level) is observed primarily within the frequency bands 0.9to 1.2 kHz 
; used in 11.4, 5.7, 3.8 and 1.9 L/h (3.0, 1.5, 1.0 and 0.5 gal/h) are shown in Figures 6.3(a) through 
ISOT 


and 2.0 to 4.0 kHz. It should be noted that within both of these frequency bands Y is not 
statistically significant at all Fourier frequencies. Figure 6.7 shows the coherence amplitude, 
Y (f), and coherence phase,@(f), for acoustic leak signals received by sensors bracketing a 


6.3(d). The hole diameters and line pressures used to establish the flow rates were 0.7 mm at 
‚the 103 kPa (15 psi) (i.e., 11.4 L/h [3.0 gal/h]); 0.5 mm at 103 kPa (15 psi) (i.e., 5.7 L/h) 
(1.5 gal/h]); 0.4 mm at 76 kPa (11 psi) (i.e., 3.8 L/h [1.0 gal/h]); and 0.4 mm at 34 kPa (5 psi) 


27 سک تا 26 


کس سا 






PERIOD (SECONDS) 


u "Wë ei 10". w* 16 








(0 f poti | WU: ፤ IM ۱ (ne xe pnmo d | WAN Fg purrr Pg 
























4 
10 (A) 3 GPH (B) 1.5 GPH 
3 
10 
z 10 
0 ۱ 1 
1 
10 pum wha 
10° AEREA ESA AR Y | 
INN ‘Il رسب‎ EC s ll | LLN CG 
Tht? f MI TILLE T MARK i 0٤١ ١ f HH tT | MITT | 32 | ] ፣ 
4 
10 (C) 1.0 GPH (ሀ) 0.5 GPH 
3 
10 
i 2 
ግሽ 
۱ 
1 
10 | 
۱ | 
10° — — — LAWANG i m a a Sud à 
سب بد اسب .سس بب‎ cru LNW لاد‎ | LL رال لل‎ 
1 2 1 2 
10 |, 6 10° 10 10 10° 
FREQUENCY (Hz) 


igure 6.3. Signal-to-noise ratio (SNR) for pipeline leaks into a sand backfill at flow rates of 11.4 L/h 
LD gah) (A), 5.7 L/h (1.5 gal/h) (B), 3.8 L/h (1.0 gal/h) (C), and 1.9 L/h (0.5 gal/h) (D). Dashed line indicates 
NR=1. SNR estimates are computed by averaging the received power at each of the sensor locations shown in 


igure 6.1. 


SNR 


Figure 6.4. Signal-to-noise ratio (SNR) for Pipeline leak into a gravel backfill. Flow rate is 11.4 L/h 
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Figure 6.5. Strength of acoustic leak signal as a function of static line pressure for a fixed hole diameter 
(0.5 mm [0.02 in]). Error bars indicate the standard deviation of six measurements used to compute the SNR at 
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Figure 6.6. Coherence amplitude and coherence phase cs a function of frequency for acoustic leak signals 
bracketing a 5.7-L/h (1.5-gal/h) leak. Sensor separation is 38.1 m (125 ft). The coherence function represents 
an ensemble average of 15 overlapping, 1024-point data segments. Dashed lines indicate 95% and 99% levels 
of statistical significance. 
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Figure 6.7. Coherence amplitude and coherence phase for acoustic leak signals bracketing a 1.9-L/h (0.5-gal/h) 
leak. Line pressure is 34 kPa (5 psi); sensor separation is 30 m (100 ft). Dashed lines indicate 95% and 99% 
levels of statistical significance. 
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coherence phase, AT), between 3.8 and 4.0 kHz for sensor pair B-C of Figure 5.1. 


Figure 6.8. Unwrapped 
line through actual data points is included. The flow rate is 11.4 L/h (3 gal/h). 
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SECTION 7 
LOCATION RESULTS 


stimate the leak location and wave speed for a given set of time 
ute the coherence function between the three sensor pairs (i.e., 

tify frequency bands of at least 100-Hz width for which the coher- 
5% level of statistical significance, (3) unwrap the coherence phase 
(4) compute the linear regression lines through each of the three 

. (4.5) through (4.7) using the known sensor positions and the com- 


wrapped phase differences between sensor pairs A-B, A-C, and 

0 to 2.5 kHz. The flow rate used in this experiment was 11.4 L/h 

t are least-squares regression ¡ines through the actual data points 
pair. The criterion for the inclusion of a phase measurement in the 
1 wave speed is that the coherence amplitude exceed the 95% level 
ach of the three sensor pairs at a given Fourier frequency. 

ilts of leak location and wave speed estimates for flow rates of 

, and 1.0 gal/h). Leak location estimates are reported as a differ- 
id actual location. 


Figure 7.1 can be used to calculate the time delays between sig- 
or pairs. The measured 00/07 values of -0.47°/Hz (A-B), 2.07°/Hz 
respond to time delays of -1.3 ms, 5.8 ms, and 7.1 ms, 


f extracting the time delays from the time series is to apply the 
Figure 7.2 shows the normalized cross-correlation coefficient, 

: (X) between time series recorded by sensors B and C. The time 

n order to isolate the high SNR portion of the leak signal spectrum 

g the correlation coefficient. Without the detailed knowledge of 

nergy provided by the coherence function, correlation analysis 

rement of the time delay between leak signals received by sensors 

correlation coefficient computed between B and C time series in 

tered from 2.0 to 2.5 kHz. Within the high-coherence interval 
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Figure 7.1. Unwrapped coherence phase between 2.0 and 2.5 kHz for sensor pairs A-B, A-C, and B-C of 
Figure 5.1. Least-squares regression lines through actual data points are included. The flow rate is 11.4 L/h 


(3 gal/h). 
used to generate th: phase curves of Figure 7.1, correlation analysis and coherence function anal- 
ysis result in approximately equal estimates of the time delay. Although this result suggests that 
the two techniques for measuring time delays are equivalent, accurate correlation analysis 
requires a priori knowledge of the frequency bands within which the acoustic leak signal is 
strong and composed of linearly propagating waves. Coherence function analysis identifies fre- 
quency bands for which the SNR is high (through the coherence amplitude) and for which the 
phase behavior is appropriate for leak location (through the coherence phase). 





Table 7.1. Leak Location and Propagation Speed Measurements 













Flow D P AP Mean Std. Mean Std. V ev Ne 
Rate’ Error Dev. Error Dev. 
(AB) (AB) (AC) (AC) 
(L/h) (mm) (kPa) (Hz) (cm) (cm) (cm) (cm) (m/s) (m/s) 
114 07 138 2100 - 2400 8.6 16.4 -2.4 23.7 1048 37 25 
114 0.7 138 3800 - 4050 18.7 29.9 14.2 31.8 917 89 18 
5.7 0.5 138. 2100 - 2400 14.4 15.8 15.8 14.9 930 136 23 
5.7 05 138 3800 - 4050 -5.8 19.8 -12.2 20.4 775 81 15 
3.8 0.5 76 3800 - 4050 -2.5 47.9 -20.7 28.1 715 150 8 











1) 11.4 L/h = 3.0 gal/h; 5.7 L/h = 1.5 gal/h; and 3.8 L/h = 1.0 gal/h 
2) 128 kPa = 20 psi and 76 kPa = 11 psi 

3) Location Algorithm Analysis Band 

4) A-B used as bracketing sensors 

5) A-C used as bracketing sensors 

6) Number of independent location estimates 
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Figure 7.2. Normalized cross-correlation coefficient, Delt, as a function of time delay (T) between time series 
recorded by sensors B and C. The time series were bandpass-filtered between 1.0 and 4.0 kHz prior to computing 
.ہم‎ The flow rate is 11.4 L/h (3 gal/h). Tac represents the predicted B-C time delay at V = 1000 m/s. 
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Figure 7.3. Normalized cross-correlation coefficient as a function of time delay between time series recorded by 


sensors B and C. The time series were bandpass-filtered between 2.0 and 2.5 kHz prior to computing p,,. The flow 
rate is 11.4 L/h (3 gal/h). Tac and Te represent predicted time delays for primary and reflected acoustic waves prop- 


agating at V = 1000 m/s. 



































































SECTION 8 
LEAK SIGNAL PROPAGATION 


The analysis of acoustic data from pipelines is complicated by the presence of n 
and multi-mode wave propagation. Multi-path signals are produced by reflections wit 
complex pipeline geometry or by signal leakage, across the connecting arms, from one 
branch of the pipeline to the other (see Figure 5.1). Multi-mode wave propagation res 
the excitation, by the leak flow field, of wave motion in different materials (e.g., gasol 
steel), or of waves in the same material that propagate at different speeds (e.g., longitu 
transverse waves). While the analysis presented above suggests that the acoustic leak 
dominated by a single propagation mode that traverses a single path from leak to 76 
imental data and simple simulations show that the effects of multi-path and multi-mod 
tion are detectable. 


The reflective nature of the pipeline is illustrated by the cross-correlation plot sh 
Figure 7.3. The primary p,, peak, which occurs at the lag time t ^ 7 ms, corresponds t 
propagating in the direction from sensor B to sensor C at speed ር ~ 1000 m/s. A secon 
which occurs at the lag time t ~ -7 ms, is consistent with reflection signals propagatir.g 
same speed, but in the opposite direction. 


Energy propagation along the pipeline results from the excitation of three types ( 
motion by the leak flow field: (1) transverse waves propagating in the steel, (2) longitı 
waves propagating in the steel, and (3) longitudinal waves propagating within the prod 
tained in the pipeline. The nominal propagation speeds for each type of wave motion ٤ 
6,000 m/s (longitudinal, steel), 3,000 m/s (transverse, steel), and 1,200 m/s (longitudin: 
gasoline). The similarity between the measured wave speed (- 1000 m/s) and the spee 
acoustic waves in gasoline suggests that in the far field of the leak, the sensors respond 
to longitudinal waves propagating through the product. These longitudinal waves are s 
indirectly through stresses induced in the steel in response to the fluctuating pressure fi 
the pipe. If other forms of wave motion are produced by the leak and are detectable, th 
measurements, and thus the location estimates, will be degraded. The detectability of 1 
nal waves propagating in steel was investigated through a calibration test in which an iı 
signal was generated by the breaking of a pencil lead near the leak location. Figure 8.1 
the time series of the calibration impulse received by sensors B and C. The measured t 
(1.2 ms) and sensor separation (7.5 m) yield a propagation speed of 6,250 m/s for the le 
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odic oscillation that occurs at approximately 50-Hz intervals with an average amplitude of 40°. 
If the signal received at each sensor is represented as a summation of a direct-path signal 
propagating at the observed wave speed and contaminating signals caused by multi-path and 
multi-mode propagation, an estimate can be made of the fraction of total energy received via the 
contaminating signals. A simple simulation in which approximately 15% of the total received 
energy was propagated by multi-path and multi-mode waves produced residual phase shifts com- 
parable to those observed in the data. 




































by the leak and are detectable in the far field, and (2) the coupling between acoustic waves in the 
product and stresses induced in the surrounding pipeline is a function of the product contained 
within the pipeline. Liquid leaks appear to be sensed primarily through energetic, low-velocity 
acoustic waves, while gas leaks are sensed via less energetic, high-velocity transverse waves 
propagating in the steel. 


edge of the impulse. This speed is consistent with the nominal value of 6,000 m/s for longitudi- 
nal waves propagated within the steel. While the calibration data do not indicate the degree to 
which the longitudinal wave mode in steel is excited by the leak flow field, it does show that 
such waves, if excited by the leak, will be detected by sensors mounted externally on the pipeline 
wall. 


SECTION 9 
PHASE UNWRAPPING 


Accurate source location requires that the location algorithm distinguish between the infor- 
mation provided by the leak signal, and ambient or system noise. The continuous nature of the 
acoustic leak signal further requires that the separation of signal from noise take place in the 
frequency domain, through coherence function analysis, rather than in the time domain. It has 
been demonstrated that source location through cross-correlation analysis is not accurate when 
applied to wide frequency bands (e.g., the 1- to 4-kHz frequency band used in Figure 7.2). 

While the location estimates given in Table 7.1 are based upon the successful application of 
coherence function analysis to relatively narrow frequency bands (100 to 500 Hz), the possibility 
exists that a similar location algorithm may be applied to frequency bands of arbitrary wi 

dth. 
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Information concerning the leak location is contained in the coherence phase. Because the 
measured phase differs from the actual phase shift by an unknown multiple of 360°, the deriv- 
ative of the coherence phase with respect to frequency, d Q/df, is required in order to estimate the 
relative arrival time of leak signals at spatially separated sensor locations. When the relative 
separation between a pair of sensors is large compared to the wavelength of the received signals, 
some form of phase-unwrapping algorithm must be applied in order to measure d ቀቆ/ቋ over a 
wide range of frequencies. Such unwrapping algorithms are easily implemented, provided that 
the coherence phase is reliably measured (i.e., the coherence amplitude is high) at many frequen- 
cies within the desired band. As the distribution of reliable phase estimates within a frequency 
band becomes more sparse, the ability to simply unwrap the phase is diminished, and the 
information provided by the phase measurements must be discarded. If the leak location and 
propagation speed of acoustic waves are known, the correspondence between measured and pre- 
dicted phase shifts can be viewed over an arbitrarily wide frequency band. 
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Figure 8.1. Time series of impulsive calibration signals recorded by sensors B and C of Figure 5.1. The estimated 
propagation speed (6250 m/s) is consistent with the nominal speed of sound in steel. 
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The excitation of transverse waves by the leak flow field, and their detectability, were 
investigated through a series of experiments in which CO,, rather than gasoline, was used as the 
product. Figure 8.2 shows the unwrapped phase shif between sensors B and C measured in the 
presence of a 103-kPa (15-psi), CO, leak. The hole diameter was 7 mm (0.029 in.). Application 
of Eq. (4.7) to the phase plot yields a propagation speed of approximately 2,400 m/s. Two 
important observations should be noted regarding this experiment: (1) the measured wave speed 
is similar to the nominal value for transverse waves propagating freely in steel, and (2) the mea- 
sured wave speed is much higher than the speed of acoustic waves propagated in CO, 

(c ~ 270 m/s). The SNR of the CO, leak was approximately 15 dB less than the SNR recorded in 
the presence of a gasoline leak at the same line pressure and hole diameter. Two conclusions 
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Figure 8.3. Unwrapped coherence phase between 2.1 and 2.4 kHz for sensor pair B-C in which the linear trend has 
been removed. The flow rate is 11.4 L/h (3 gal/h). 
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Figures 9.1 and 9.2 show the unwrapped phase shift between sensors A-B, A-C, and B-C 
in which the unknown multiples of 360° required to unwrap the phase were computed from the 
predicted $(f) lines (shown as solid lines in the figure). Reliable phase measurements (indicated 
by markers in the plots) correspond to coherence amplitudes that exceed the 95% level of statis- 
tical significance; the flow rates are 11.4 L/h (3 gal/h) (Figure 9.1) and 5.7 L/h (1.5 gal/h) 
(Figure 9.2). The frequency distribution of reliable phase measurements for the 11.4-L/h 
(3-gal/h) data is such that all of the information contained in the 2- to 4-kHz band can be used in 

















Figure 8.2. Unwrapped coherence phase between 2.7 and 3.0 kHz for sensor pair B-C of Figure 5.! in which CO, is 
used as the product. The line pressure is 103 kPa (15 psi); the hole diameter is 0.7 mm. The estimated propagation 
speed (2400 m/s) is consistent with the nominal speed of transverse waves in steel. 

















The effect of multi-path and multi-mode wave propagation can also be observed in the 
coherence phase. Figure 8.3 shows a plot of the phase shift between sensors B and C in which 
the linear trend has been removed. The residual phase shift is dominated by a non-random, peri- 






may be drawn from these measurements: (1) freely propagating transverse waves are produced 
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PHASE SHIFT (DEGREES) 





location estimate ከሃ implementing a straightforward phase-unwrapping algorithm. As the 
ሃ rate 15 reduced, however, the simple unwrapping algorithm works only within a small num- 
of narrow frequency bands (e.g., 2.2 to 2.5 kHz, and 3.7 to 4.0 kHz in Figure 9.2). The simi- 
y between the measured and predicted phase shift outside of these narrow bands suggests 


a more robust unwrapping algorithm may be capable of exploiting a greater fraction of the 
lable phase information for the purpose of leak location. 
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e 9.1. Unwrapped coherence phase between 1.5 and 4.5 kHz for sensor pairs A-B, A-C, and B-C. Solid lines 
عا‎ predicted coherence phase for linearly propagating plane waves based upon known leak location and propa- 
speed. Flow rate is 11.4 L/h (3 gal/h). 
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Figure 9.2. Unwrapped coherence between 1.5 and 4.5 kHz for sensor pairs A-B, A-C, and B-C. Solid lines 
indicate predicted coherence phase for linearly propagating plane waves based upon known leak location and propa- 
gation speed. Flow rate is 5.7 L/h (1.5 gal/h). 
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